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We generalize the Gutzwiller approximation scheme to the calculation of nontrivial matrix el- 
ements between the ground state and excited states. In our scheme, the normalization of the 
Gutzwiller wave function relative to a partially projected wave function with a single non projected 
site (the reservoir site) plays a key role. For the Gutzwiller projected Fermi sea, we evaluate the 
relative normalization both analytically and by variational Monte-Carlo (VMC). We also report 
VMC results for projected superconducting states that show novel oscillations in the hole density 
near the reservoir site. 



I. INTRODUCTION 

This paper concerns the calculation of matrix elements 
using projected wave functions of the form, Y$>) = P . 
Here, P — Yii(^ ~ n i^ n il) is a projection operator which 
excludes double occupancies at sites i, and l^o), a trial 
wave function. Projected wave functions of this form 
were originally proposed by Gutzwiller to study elec- 
tronic systems with repulsive on-site interactions 1 . The 
choice of |^o) depends on the problem under considera- 
tion. For instance, a projected Fermi liquid state, 

^i*fs) = p n c Mj°> • w 

was used successfully in the description of liquid 3 He 
as an almost localized Fermi liquid 2,3 . Soon after the 
discovery of high temperature superconductivity in the 
cuprates, projected BCS wave functions were proposed 
as possible ground states of the so-called t — J models. 
Early results from variational Monte Carlo (VMC) stud- 
ies as well as a renormalized mean field theory based 
on Gutzwiller approximation showed that a projected d- 
wave BCS state, 

PnP |*bcs) = PnP J[ (u k + u*4 T c -ki) 1°) > ( 2 ) 
k 

reproduces many features seen in the phase diagram of 
the high temperature superconductors&LiV The projec- 
tion operator FV which fixes the particle number N in 
is useful when considering the phase diagram near 
half filling 6 .. Without Pn in @, one would need to con- 
sider the effects of particle number fluctuations, which 
become singular near half-filling 910 . 

Detailed VMC studies have been carried out recently 
using projected d-wave BCS states as variational wave 
functions for the two dimensional Hubbard modelii, after 
a suitable canonical transformation^. Similar wave func- 
tions have been proposed in the literature for cobaltate 
superconductors as well as organic superconductors 1 ^* 1 *?-. 
To make analytical progress however, it is desirable to ex- 
tend Gutzwiller 's scheme and construct normalized sin- 
gle particle excitations and calculate matrix elements. In 
this paper, we take the first step in this direction. Wc 



construct normalized excitations of the Gutzwiller pro- 
jected Fermi sea and consider the evaluation of matrix 
elements. 

In his original paper, Gutzwiller proposed that in cal- 
culating expectation values of operators with projected 
wave functions, the effects of projection on the state |^o) 
could be approximated by a classical statistical weight 
factor, which multiplies the quantum resullpii Thus, for 
example, 

(gjOjg) ^ (*o|Q|*o) ^ 
~ 9 (* |*o> ' 1 ' 

where O is any operator, and g, a statistical weight fac- 
tor. The basic idea is that the projection operator P re- 
duces the number of allowed states in the Hilbert space, 
and invoking a simple approximation, such a reduction 
can be taken into account through combinatorial factors. 
For example, expectation values of the the kinetic energy 
operator c|cj+cjcj and the superexchange interaction be- 
tween sites i and j, Si ■ Sj in the projected subspace of 
states are renormalized by the Gutzwiller factors, 

9t ~ l-n/2' 9s = (l-n/2) 2 ' (4) 

where n is the density of electrons. In deriving these 
renormalization factors, one considers the number of 
states that contribute to (*|0|*) and to (W |6|* ) re- 
spectively. The ratio of these two contributions is iden- 
tified as the renormalization factor. 

It is clear that this approach can be generalized to 
evaluate matrix elements of an operator O between dif- 
ferent projected states. However, as we will see in this 
paper, many of the matrix elements that are of interest 
are reduced to the calculation of matrix elements between 
partially projected wave functions of the form 

= p( i*o>, p/=na- • (5) 

The wave function \^',) describes a state where double 
occupancies are projected out on all sites except the site 
I, which we call the reservoir site. The reason for the ap- 
pearance of reservoir sites is not far to seek. Consider, for 
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hole doping 



Double occupancy is projected out on all sites except the 
site I, called the reservoir site. Unless specified other- 
wise, we take |\&o) to mean the Fermi sea. For the calcu- 
lation of single particle excitations and matrix elements, 
we need expectation values such as 

(nd\*'i) = , (tto|6|*o) m 

WW) 9 <* |*o> ' 1 ' 

that generalize the Gutzwiller renormalization scheme 
© to partially projected wave functions. 



FIG. 1: Double occupancy of the reservoir site, (^t^U)*' — 
1 — X, as a function of doping, for the partially projected 
Fermi sea. See @ and Q . Note the good agreement between 
the Gutzwiller result (solid line) , Eqs. 111)1151 and the VMC 
results for the projected Fermi sea (open circles). Statistical 
errors and finite-size corrections are estimated to be smaller 
than the symbols. 



example, the operator Pcjf. Clearly it can be rewritten 
as c\\P[. Since calculation of matrix elements involv- 
ing excited states involve the commutation of projection 
operators with creation/destruction operators, partially 
projected states arise inevitably within the Gutzwiller 
scheme. 

In this paper, we present a method to calculate matrix 
elements between a partially projected Fermi sea, i.e., 
a projected Fermi sea with a reservoir site at I, as in 
(J5J. We will show that this problem has to be solved if 
we were to construct normalized particle/hole excitations 
of the (fully) projected Fermi sea. The same problem 
arises when calculating matrix elements for particle/hole 
tunneling into the projected Fermi sea. We develop an 
analytical approximation to solve this problem, and use 
it to calculate various matrix elements. We use VMC to 
test the validity of the approximation and find that our 
analytical results for the partially projected Fermi sea 
are in good agreement with the results from VMC. 

The outline of the paper is as follows. In Sec. [HJ we 
present results for the occupancy of the reservoir site. We 
use these results in Sec. IIIII where we show how normal- 
ized single particle excitations can be constructed from 
the projected Fermi sea. In Sec. IIVI we calculate the 
matrix elements for particle/hole tunneling into the pro- 
jected Fermi sea. VMC results for density oscillations in 
the vicinity of the reservoir site for both projected Fermi 
sea and BCS states are presented in Sec. |V] The final 
section contains a summary and discussion of results. 



II. OCCUPANCY OF THE RESERVOIR SITE 



Consider a partially projected wave function, 

= P/|* ), P{ = - n iT ni;) . (6) 



A. Gutzwiller approximation 

In order to evaluate the generalized renormalization 
parameters <?' in J7J), we obviously need the normalization 
{%\%). We define 



X 



(*o\PP\*o) 
(*o|P/P/|*o) 



mm) 



(8) 



the norm of the fully projected state relative to the state 
with one reservoir site. Invoking the Gutzwiller approxi- 
mation, we estimate this ratio by considering the relative 
sizes of the Hilbert spaces, 
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L! 

N r W±lN h \ 



LI 



(9) 



JVf ! AT j ! Nh ! t (JV T -l)!(JVi-l)!(iV h + l)! 



where L — iVj + iVj, + Nh, is the number of lattice sites, 
Vj, N± and Nh, the number of up spins, down spins and 
empty sites respectively. The first term in the denom- 
inator of © represents the number of states with the 
reservoir site being empty or singly occupied; the second 
term represents the state with the reservoir site being 
doubly occupied. 

Eq. can be simplified in the thermodynamic limit. 
We get, 



X 



1 



(l-n T )(l-n;) ' 



(10) 



where the particle densities, n a = N a /L (a = f,|) and 
n = n-f +«x- The above argument can be extended to the 
case of two unprojected sites in an otherwise projected 
Fermi sea. We then get, 



(* |PP|*0) = X 2 



(*o|P/ m PLI*o) 



(11) 



where, P lm = Yli^mi 1 ~ n «,T n u)- We note for later use 
that 



1 — X ntn 



1 n l 



X (1 - n) 



(12) 
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B. Exact relations 

Assuming translation invariance, it is possible to derive 
the following exact expressions 

((l-n, T )(l-n u ))*j - X(l-n) (13) 
(ni a (l - n;_ CT ))*; = In, (14) 
(<£)*; = (n; T n u )*j = 1 - AT (15) 

for the occupancy of the reservoir site, where 

(„.),{ = (^...m/in^'i) ■ 

The proof is straightforward. Consider for instance, the 
probability l|13|l of finding the reservoir site empty. Since, 

(9 \P(l-ni)P\^ ) (16) 
= (* |P/(l-n iT )(l-n u )P/|*o) 

we have, 

((1 -n n ){l -niy))^ 

Eqs. I|14|l and l|15[l can be proved analogously. 

C. VMC results for projected Fermi sea and BCS 
states 

In Fig. ^ we compare (|10|) with VMC results for 
(d)^> = 1 — A. We find that the results from the general- 
ized Gutzwiller approximation are in excellent qualitative 
agreement with the VMC results for a partially projected 
Fermi sea. We also used VMC to obtain the same quan- 
tity using projected s/d-w&ve BCS states as variational 
states in the simulation. The results for (d)^> in BCS 
states are shown in Fig. |2 In contrast to the projected 
Fermi clear deviation from the Gutzwiller approxi- 

mation is seen. This underscores the importance of pair- 
ing correlations in the unprojected wave function that 
are not completely taken into account by the Gutzwiller 
approximation scheme. These differences between Fermi 
sea and BCS states are discussed in more detail in Sect. 
Ivl where we consider density oscillations in the vicinity 
of the reservoir site. 

In the following we discuss some details of the VMC 
calculations with one unprojected (reservoir) site I. As 
mentioned earlier, single occupancy is enforced (by pro- 
jection) on all other sites. Simulations are performed on 
a finite square lattice spanned by two vectors (L x ,L y ) 
and (—L y , L x ) with periodic boundary conditions^. The 
number of sites, L = L\ + L y . The numbers of up- and 
down-electrons are chosen to be equal, iV-j- = iVj_. The 
simulation for the local quantity (d)-q,> — (ni^riii)^' has a 
larger statistical error than results for macroscopic quan- 
tities in uniform systems because the summation over site 



indices yields effectively L times more statistics for the 
latter. In order to overcome this problem, we update 
the reservoir site more often than the projected sites. 
Accordingly, the transition probability needs an extra 
weighting factor to keep the local balance. With this pro- 
cedure, we can improve the statistical accuracy by about 
one order of magnitude. In addition, we carry out mea- 
surements after every update. Usually, in VMC simu- 
lations, measurement are performed every O(L) updates 
to obtain independent samples since similar states return 
similar sampled data. However, in the case of n^rtyj,, a 
measurement returns only or 1; viz., the sampled data 
can be different even when the states are similar. Given 
this, the measurement after every update seems more 
reasonable as it reduces statistical errors. Furthermore, 
we have restricted updates to the transfer of a single elec- 
tron to an unoccupied site, and excluded updates via the 
exchange of two electrons. The calculation of the tran- 
sition probabilities for the former update consumes time 
of 0(N a ) whereas the time taken for the latter update 
is 0(N%). As the system size increases, this restriction 
achieves efficiency. We have collected statistics from up 
to 60 independent runs over two days, and the total num- 
ber of updates amounts to 10 8 ~ 10 9 . 

For superconducting states, one can perform the VMC 
simulation either with fixed particle number PnP\^bcs) , 
or with a fixed phase P|\E , BCs)^ 0,i ^- For the latter choice, 
particle number fluctuation hinders the variational wave 
function from reaching half filling unless the chemical po- 
tential fi goes to infinity. On the other hand, the wave 
function can be optimized by varying the gap even 
at half filling, if we choose to fix the particle number. 
It is important to note that simulations with fixed par- 
ticle number are done not with the most probable N of 
PI^bcs), but that of I^bcs)- This is because P decreases 
the average particle number 10 . Despite these differences, 
both choices of wave functions yield quantitatively sim- 
ilar results^. Throughout this paper we choose to fix 
the particle number while working with projected BCS 
states. 

Let us define = v^/uk- For the d-wave BCS state, 
ak—o — in the thermodynamic limit. However, if one 
chooses a,k=o = in the finite system, the fc = state is 
unoccupied, although it is the lowest energy state. One 
can also choose a large value for ak=o- Usually, the dif- 
ference between these choices is 0(1/N). We expect an 
N electron system with with ak=a = to be similar to 
an N + 2 electron system with large a,k=o, because the 
two electrons that are no longer in the k = state can 
occupy other available states. However, this argument 
fails at half filling for projected states, because there are 
no available states left for the two extra electrons. So 
it should not be surprising that X depends strongly on 
these choices close to half filling; the dfc = o = definition 
gives larger X than the other does as shown in Fig. [21 At 
the other fillings, our results show only 0(1 /N) of differ- 
ence between these choices. Except for Fig. [21 where we 
show both cases, all other results in this paper are ob- 
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hole doping 

FIG. 2: Double occupancy at the reservoir site (nt n l)*' = 
1 — X, as a function of doping, for the partially projected 
BCS wave function, see © and (|5J. The parameterization 
follows Ref. IHi Statistical errors and finite-size corrections 
are estimated to be smaller than the symbols. 



III. SINGLE PARTICLE EXCITATIONS OF THE 
PROJECTED FERMI SEA 

We consider the particle excitation 

l*2a> = ^4.1*0) , (17) 

and the hole excitation 

l*to> = Pck^o) ■ (18) 

Any calculation involving l^^.) needs the respective 
norms, N^ a = {^^a^ka)- We now calculate these norms 
within the generalized Gutzwiller approximation. 

A. Particle excitation 



tained for a choice of large ak=o', i-s., we take aj, =0 larger 
than any other a^. 

The system size dependence is quite small except in the 
vicinity of half filling. In fact, it is qualitatively consis- 
tent with the Gutzwiller approximation; size dependence 
enters only as (Nh + 1)/L in Eq. © and is negligible 
for large Nh- In Fig. [3J we show the dependence of (d) 
on the system size. As shown in Fig. [31 (d) approaches 
unity for the projected Fermi sea. For the projected d- 
wave BCS state, we speculate that the value of (d) goes 
to unity too, because it does not saturate, but increases 
more rapidly as 1/L decreases. 



For the particle excitation, we get, 
N + 

j^f = 1 - n + gt K - n t) =9t(l- «L) , (19) 

where g t = (1 - n)/(l - n a ), N G = (#|#), and n\ a = 
( c lo- c fco-)*o i s the momentum distribution function in the 
unprojected site. 

Equation i|19f) has appeared frequently in the litera- 
ture, here, we repeat its derivation to facilitate a com- 
parison with the analogous problem for hole excitations. 
The norm (^ka-l^kcr) * s gi ven by 



= (^\c ka PPcl\^ ) = i^e ife C— )(*o|P/(l-^- CT )c iCT cL(l-n m - <T )P4l*o) 

l,m 

i^ol-P/Cl-n^Cl-n^P/l^o) + j £ e^- m H^ \Pc lrT cl a P\^ Q ) 



(m-nm N Gsr ik(l _ m) (V\cl«cM 

<*|*> L (*|*> ' [M) 

where we have used (|16H for the diagonal contribution in the last step. Invoking the Gutzwiller approximation for the 
off-diagonal term, Eq. (|19|) follows directly from l|2(JI) . 



B. Hole excitation 



The normalization of the hole excitation can be done analogously. We get, 

where, P; TO = Yii^i m Q ~ n iA n i,i)- The last term in the above equation corresponds to a hopping process between 
two reservoir sites. The generalized Gutzwiller approximation assumes that the matrix elements are proportional to 
the square roots of the corresponding densities (|13I14I15|I . 



+ wrT, elk{l ~ m) ^o\PLclc m M* ) , 
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Invoking the Gutzwiller approximation and using we get, 



N, 



ha 



N G 



1-X 
X 



X 2 n a (l - n a ) 



^X{l-n)^Xn~+^fXn~ZVT^~X 



(21) 



for the normalization of the hole excitation relative to the norm of the Gutzwiller wave function. 

The general expression ()21(l for the hole normalization, can be simplified upon using the Gutzwiller result ljl2|l for 
the relative norm X . We then get, 



VT^^+^T^^{i-x)/x 



n a (l - n a ) 

for the last term in l|21|) . Finally, we get the simple result, 



(n 



o 

ka 



[(1 - n) + n_ a ] 2 
(l-n CT )(l-n) 



1 - n-g 



ly k,a 1 - n- a n ka 

-T7— = nka-p. r = 22) 

N G (1 - n) g t 

It is interesting to compare this result for the normaliza- 
tion of the hole excitation with the corresponding expres- 
sion (|19|) for the particle excitation. The vanishing of the 
latter at half filling could have been expected. But the 
divergence of N^ a as n — > 1 is surprising. We will return 
to this point in the next section. 



IV. TUNNELING MATRIX ELEMENTS 

We now consider the tunneling of electrons and holes 
into a projected wave function. Single particle tunneling 
into a projected superconducting state has been consid- 
ered recently by Anderson and Ong 9 , and Randeria et 
ali^. Here, we restrict ourselves to the projected Fermi 
liquid state and evaluate the tunneling matrix elements 
by retaining systematically, all terms arising from the 
commutation of the electron creation and destruction op- 
erators with the projection operator P, as outlined in Sec. 



C. Consistency check 



The norm Nj^, has to vanish whenever |\& 



ko 



Pc ktT \^>o) vanishes. For the Fermi sea this is the case 
when k < 1zf, i.e. when vP ka = 1. This physical condition 
is obviously fulfilled by 119fl . Similarly, we expect to 
vanish for n% a — 0, which is satisfied by l122[| . Thus, the 
Gutzwiller result llUfl obeys the normalization condition 
for the hole excitation and the theory is consistent. 
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FIG. 3: The size dependence of the double occupancy at the 
reservoir site, 1151 . at half filling. The VMC result for the 
Gutzwiller state (upper curve) seems to converge nicely to 
unity in the thermodynamic limit, in agreement with Eqs. 
HUH and (11 -it . The results for the projected d-wave show a 
pronounced dependence on the occupancy of the k — state. 
See text for details. 



A. Particle tunneling 

Consider first, the matrix clement 



ML = 



+ \ j 



LI*) 



(23) 



The numerator may be calculated easily by using the 
result of (JTSJl: 



*o|c fcg P C L-P|*o) 
N G 



N G 

= st(i-nL) • 

From the above expression we find that the particle 
tunneling matrix element takes the form, 

It vanishes at half filling n — > 1, implying that the ad- 
dition of electrons is not possible exactly at half filling 
because of the restriction in the Hilbert space. 

B. Hole tunneling 

Next we evaluate the matrix element 

2 



N^Ng 



(25) 
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corresponding to the tunneling of holes into the projected 
state. Naively, we might expect this process to be allowed 
at half filling, since the removal of electrons is not for- 
bidden by the projection operator. Consider now, the 



matrix element in the numerator of 125(1 . We follow the 
same procedure used to evaluate the norm of the hole 
wave function in Sec. IIII Bl and use ((12(1 and 1(14(1 and 
find, 



<*o|cL-Pc fc(T P|*o) 



N G 
Xn 



— E< 



x ik{l—rn) i 
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Jr. J, ^ 
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l,m l^m 
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[^/x^ZVT^x + y/X(l - n)^Xn^] vT 



X(l - n a )n a 



(26) 



= n a + (n° ka - n a ) 



ii-aUg + (1 - n)n a 
(1 - n a )n a 



'Act 



Using this expression together with the norm 1(21(1 of the 
hole excitation, we obtain the hole tunneling matrix ele- 
ment P5)l. 



= n irrr = (27) 



a surprising result, in that it vanishes at half filling 
(n-f = n-i = 0.5) too. 

The vanishing of the hole tunneling matrix element at 
half filling is clearly related to the divergence of the norm 
of the hole excitation. This, in turn, is related to the fact 
that X — > 0, as n — > 1 (cf. Ea lfTU|) '). The vanishing hole 
tunneling matrix element can then be understood as fol- 
lows. When the reservoir site is doubly occupied, a single 
hole in the otherwise projected Fermi sea can be found 
in any of the lattice sites. Consequently, when double 
occupancy of the reservoir site occurs with probability 1 , 
as it does at half filling, an "orthogonality catastrophe" 
occurs leading to zero overlap for the tunneling matrix 
element. Note that the result ((27(1 hinges on the exact 
functional dependence l(12|) of (1 — X)/ X on the particle 
densities n a . On the other hand, the particle tunnel- 
ing matrix element Mi is not affected by the functional 
form of the relative normalization factor X. If X were to 
vanish more slowly than (1 — n) at half filling, then from 
((21(1 and ((26(1 , one could conclude that the hole tunneling 
matrix element M^ a does not vanish as n — > 1, possibly 
leading to an asymmetry between particle and hole tun- 
neling. Our analytical results preclude this possibility 
for the projected Fermi sea. But we are unable to pro- 
vide a definite answer for the projected superconducting 
states, in view of the discrepancy between the Gutzwillcr 
approximation and the VMC results (Fig. EJ) . To under- 
stand this discrepancy, we study density oscillations in 
the vicinity of the reservoir site using VMC. 



V. DENSITY OSCILLATIONS NEAR THE 
RESERVOIR SITE 

To clarify the limitations of the Gutzwiller approxima- 
tion for projected superconducting states, we use VMC 
to calculate the hole density in the vicinity of the reser- 
voir site. We find that the density oscillations seen are 
very different for the projected Fermi sea and the BCS 
states. 

Fig. 01 shows VMC results for the hole density 

n h (m) = (1 - ra m )*; , 

in the partially projected state j'J/j) are presented in the 
first row of. The sites to are distinct from the reser- 
voir site I (marked by a cross in the figure). All results 
shown correspond to half filling; viz., — = 0.5. 
We choose A = 1 for the BCS states. The vectors 
of periodic boundary conditions are L\ = (L x ,L y ) and 
L2 = (—L y ,L x ) respectively, with L x = 37, L y — 1; In- 
cluding the reservoir site, L = L? x + l? x — 1370 sites. 
In the figure, white/black correspond to high/low val- 
ues of n^(rn), which is scaled by a logarithmic gray scale 
varying in the range —8.5 < log (to) < —6. Thus, the 
same gray represents the same value in all the three cases 
shown. 

For the Fermi sea, we see that the hole is distributed 
more uniformly than the other cases even though the 
diagonal direction has a larger probability of being oc- 
cupied by a hole. The s-wave shows a checker-board 
pattern. The ci-wave has a quasi checker-board pattern 
where only one of four sites is black, and the hole tends 
to be near the reservoir site. The VMC results for the 
projected BCS wave functions are strikingly different in 
that the hole density is not uniform. On the other hand, 
the Gutzwiller approximation would be exact, if all states 
in the Hilbert space contribute equally to the wave func- 
tion. That would correspond to a uniform density of 
holes. Clearly, the Gutzwiller approximation has to be 
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extended to treat projected superconducting wave func- 
tions. This is in agreement with our previous consider- 
ations, where we found that the functional form of X 
(Eq. 1101 derived using Gutzwillcr approximation) agrees 
with the VMC calculations only for the projected Fermi 
sea, but not for BCS states (see Fig.^and FigEJ. 

To further investigate the effects of Gutzwiller projec- 
tion, we also plot (second row of Fig. the correlation 
function 

d h\ m ) = { n l-\nn{l - n mT )(l - n m ;)) (28) 
- {ni\n n ){(l-n m] )(\-n mi )) 

in systems without the Gutzwiller projection. This cor- 
relations function between a hole at site to and a doubly 
occupied site at I corresponds to the quantity nh(m) for 
the partially projected wave function close to half filling. 
This is because, in the latter case, the unprojected site 
is doubly occupied. Note that translation invariance im- 
plies that the second term in (|28|l does not depend on 
the site indices / and m, and is a constant factor. Then, 
using Wick's theorem, the correlation function d ® (to) 
is reduced to a function of (c\^Cj^) and (ci,tCj,l)- The 

quantity di ' (to) can be evaluated exactly, performing 
a Fourier transform (we use the same system size and 
boundary conditions). The logarithm of the correlation 
function is scaled in the second row of Fig.0|by the gray 
scale varying in the range (—22, —4). Both (cj jCj.f ) and 
( c i,1 c j,l) show Friedel oscillations. 

For the Fermi sea, only (c^Cj^) is finite. The nesting 
of the Fermi surface by Q — (tt,tt) then leads to the to 
the checker-board pattern for the hole density observed 
in Fig. 0] (second row). 

For the s-wave, the Friedel oscillations of (c\^Cj^) are 
similar to that of the Fermi sea while the oscillation of 
( c i,1 c j-l) 1S pbase shifted by ir/2. Summing both contri- 
butions to d ® (to) the oscillations are smeared out. In 

contrast, for the <i-wave both (cj *Cj,t) and (ci^Cj^) os- 
cillate in phase, leading to the oscillation observed. 

Let us compare these results with those obtained after 
projection. For the Fermi sea, we see clearly that the 
density oscillations are suppressed by projection. This is 
likely because projection reduces the discontinuity at the 
Fermi level, thereby suppressing the nesting by Q and 
the corresponding Friedel oscillations. 

The emergence of the checker-board pattern in the pro- 
jected s-wave suggests that Gutzwiller projection affects 
( c i,t c 3,l) stronger than (c\^Cj^). With only one contri- 
bution, the Friedel oscillations are no longer smeared out 
and are observed. 

Projection changes the pattern qualitatively for the d- 
wave too. The observed pattern resembles approximately 
the function, ~ sin 2 (a;7r/2) sin 2 (?/7r/2) (see Fig. 0J top 
row), with to = (x,y). This indicates that the nodal 
points at (±^,±7^) contribute dominantly after projec- 
tion. Furthermore, in this case, the hole tends to stay 
near the reservoir site. It means that only a part of the 



Hilbert space has a large weight, leading to a deviation 
from the Gutzwiller approximation. We believe this ef- 
fect cannot be captured within the Gutzwiller approxi- 
mation without invoking off-site correlations^. 



VI. DISCUSSION 

In this paper, we extended the Gutzwiller approxi- 
mation scheme to construct normalized excitations and 
matrix elements for the projected Fermi sea. In typi- 
cal calculations, one needs to determine matrix elements 
between partially projected Gutzwiller projected states, 
where double occupancies are projected out at all but one 
site I (called the "reservoir" site) . The occupancy of the 
reservoir site, n\ turns out to be an important quantity in 
the calculation of matrix elements. Since the wave func- 
tion projects out double occupancies on all sites to 7^ I, it 
follows that the occupancy n m < 1. But, ni = {0, 1,2}. 
Therefore, our results for ni are nontrivial in that the 
Gutzwiller approximation is extended to calculate the oc- 
cupancy at an unprojected site. We presented an analyti- 
cal method to calculate such matrix elements and showed 
that the approximations are in good agreement with re- 
sults from variational Monte Carlo (VMC) for the Fermi 
sea. These results were used to construct normalized sin- 
gle particle excitations of the fully projected Fermi sea, 
and to calculate matrix elements for tunneling into the 
projected Fermi sea. 

Single particle tunneling in projected BCS wave func- 
tions has been discussed recently, by Anderson and Ong 9 , 
and Randeria et ali^. In our calculations for tunneling 
into the projected Fermi sea, we find the surprising result 
that the matrix elements for both particle and hole tun- 
neling vanish as n — > 1 (half filling) . Within our scheme, 
the result follows from the behavior of the charge den- 
sity in the vicinity of the reservoir site. In particular, for 
the projected Fermi sea the analytical result hinges on 
the expression for single occupancy of the reservoir site, 
Eq. (|10fl . As can be seen in Fig. 2, the analytical result 
does not agree with numerical calculations done for pro- 
jected BCS wave functions. This discrepancy underscores 
the importance of pairing correlations in the unprojected 
wave functions, which are not taken into account within 
the Gutzwillcr approximation scheme. 

There are two ways by which electron correlations arise 
in the Gutzwiller scheme: one is through the mean field 
or trial wave function l^o), an d the other via the pro- 
jection on the subspace of no double occupancy, |\E') = 
P\^fo). The latter effect, which results in the reduction 
in the size of the Hilbert space can be described by com- 
binatorial arguments, leading to (|10fl . As seen in Fig.^ 
the analytical and VMC results are in good agreement 
for the case of the projected Fermi sea. We can trace 
this agreement back to the fact that the Fermi sea does 
not contain any additional explicit correlations. 

Consider instead I^bcs)) which contains additional, 
molecular field correlations in the unprojected wave func- 
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FIG. 4: Top row: VMC results for the hole density n/,(m) = {(1 — 7i m ))if (color coding: white/black correspond to high/low 
values of nh(ra)) in the partially projected state \^{), for sites m other than the reservoir site I (marked by the cross). Left: 
Fermi sea. Middle: s-wave state. Right: d-wave state. Second row: Exact results for dh(m), see Eq. 1281 . in the unprojected 
state (color coding: white/black correspond to high/low values), for sites m other than the doubly occupied site I (marked by 
the cross). Left: Fermi sea. Middle: s-wave state. Right: d-wave state. All: The lattice has 37 2 + 1 = 1370 sites and n = 1. 
For the BCS states, A = 1. 



tion. Here, we may expect deviations for quantities like 
the relative normalization X from the combinatorial re- 
sult iJTDJl. Indeed, the VMC data presented in Fig. H 
confirms this expectation. For instance, the data show 
a qualitatively different dependence of X on doping, for 
the s-wave BCS states. The VMC data indicates a pos- 
sibly different limiting behavior for X in the limit n — > 1, 
as indicated by the analysis of the data as a function of 
inverse cluster-size, presented in Fig. |3| 

For the s-wave BCS state, we observe a dramatic en- 
hancement in the double occupancy at the reservoir site 
for low doping, which we understand as a consequence 
of enhanced on-site pairing, relative to the Fermi liquid 
state. On the other hand, The double occupancy of the 
reservoir site is reduced for the <i-wave, since the d-wave 
state suppresses on-site pairing fluctuations. The quanti- 
tative behavior of the normalization ratio X as a function 
of doping, for projected superconducting states is thus a 
subtle problem which we hope to solve in the future. 

We also studied the hole density near the reservoir site 
for projected superconducting wave functions at half fill- 



ing using VMC. The results are shown in the top row 
of Fig. 0] For the projected Fermi sea, we find that the 
hole density is uniform. However, for the superconduct- 
ing states, we find that projection induces oscillations in 
the hole density near the reservoir site. For the projected 
d-wave state, we find that the hole density is mostly near 
the reservoir site. We believe that the Gutzwiller ap- 
proximation needs to be extended to treat pairing cor- 
relations in the superconducting wave functions to un- 
derstand these results fully. This issue, along with the 
study of systems away from half filling and their possi- 
ble relevance to the checker-board pattern observed in 
scanning tunneling microscopy of the high temperature 
superconductors is left to future research. 
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